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Abstract 


We propose a simple texture of the neutrino mass matrix with one sterile neutrino 
along with the three standard ones. It gives maximal mixing angles for Ue —> us 
and —> Ut oscillations or vice versa. Thus with only four parameters, this mass 
matrix can explain the solar neutrino anomaly, atmospheric neutrino anomaly, LSND 
result and the hot dark matter of the universe, while satisfying all other Laboratory 
constraints. Depending on the choice of parameters, one can get the vacuum oscillation 
or the large angle MSW solution of the solar neutrino anomaly. 


Recently the super-Kamiokande experiment has conhrmed the atmospheric neutrino os¬ 
cillation result, suggesting nearly maximal mixing of with another species of neutrino 
The same experiment has also conhrmed the solar neutrino oscillation result, which 
suggests mixing of z/g with another species of neutrino |Q. Moreover, the energy spectrum 
of the recoil electron seems to favour the large mixing-angle vacuum oscillation of z/g over 
the MSW solutions |^, although this may have limited statistical signihcance in the global 
ht to the solar neutrino data They have led to a hurry of phenomenological models 

for neutrino mass and mixing which can account for these oscillations [5-8], most of which 
are focussed on the bi-maximal mixing angles for the atmospheric and the solar neutrinos. 
However, almost all of these works are based on the three-neutrino formalism, involving the 
standard left handed neutrinos z/g, z/^ and Vr IH 1 - 


On the other hand the inclusion of the LSND neutrino oscillation result 0 is known 
to require a fourth neutrino, which has to be a sterile one (z/ 5 ) for consistency with the 
observed Z-width [|T^]. Moreover it requires either z/^ or z/g to oscillate into z /5 for explaining 
the atmospheric and solar neutrino anomalies, while requiring z/^ —> z/g oscillation for the 
LSND result. Thus the three-neutrino models for atmospheric and solar neutrino anomalies, 
based on a z/g — z/„ — mixing, are in direct conflict with the LSND result. While the LSND 


result has not been corroborated by the preliminary KARMEN data |]^, the statistical 
signihcance of the latter is limited by its lower sensitivity in the relevant region of parameter 
space. Indeed, with the standard statistical method the 90 % c.l. limit of KARMEN excludes 
only half the parameter space of the LSND data in the Am^ < 2eV‘^ region |]T^. Hopefully, 
this issue will be resolved by the proposed mini-BOONE experiment at Fermilab along with 
more data from KARMEN. It seems to us premature, however, to rule out the LSND result 
at present. Therefore we have tried to construct a four-neutrino mass matrix, which can 
account for the present solar and atmospheric neutrino data along with the LSND result. 
It can also account for the hot dark matter content of the universe |^, while satisfying all 
laboratory and astrophysical constraints |T^, |T^, . 


Table 1 summarises the experimental constraints on neutrino mass and mixing param¬ 
eters, which are relevant for our model. The large angle MSW and the vacuum oscillation 
solutions to the solar neutrino data B O, are taken from a recent ht to the z/g sup¬ 


pression rates along with the recoil electron spectrum by Bahcall, Krastev and Smirnov 
0. For both the solutions the ht favours the oscillation of z/g into a doublet neutrino over 
^s- The reason is that in the former case the NC scattering of this doublet neutrino 
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Table 1: Present experimental constraints on nentrino masses and mixing 


Solar Neutrino 0 

Am? ~ (0.8 — 2) X 10 

(Large angle MSW) 

sin^ 20 ~ 1 

Solar Neutrino 

Am^ ~ (0.5 - 6) X 10-^°eV2 

( Vacuum oscillation) 

sin^ 20 ~ 1 

Atmospheric Neutrino |]I| 

Am^ ~ (0.5 — 6) X 10“^eV^ 
sin^ 20 > 0.82 

LSND 0 

Aml^ ~ (0.4 - 2)eV^ 
sin^ 20 g^ ~ 10-3 - 10-2 

Hot Dark Matter |13 

E.m, ~(4-5)eV 

Neutrinoless Double Beta Decay |14 

< 0.46eV 

CHOOZ p 

Amlx < 10~^eV^ 

(or sin^ 20 gx < 0 . 2 ) 


with electron can partly acconnt for the discrepancy between the observed suppression rates 
in super-Kamiokande and the Homestake experiments. On the other hand one can get ac¬ 
ceptable solutions with z/g —us oscillation if one makes allowance for a 20 % normalisation 
uncertainty for the Homestake experiment. This will also enlarge the acceptable range of 
Am^. Therefore we shall consider both the oscillation scenarios z/g —^ z /5 and z/g ^ Ur in 
our model. It should be added here that the best value of sin^ 29 for the large angle MSW 
solution is slightly less than 1 ; and even there the quality of fit is rather poor when all the 
experimental data are put together 0. However one can get acceptable £t with the large 
angle MSW solution, including the sin^ 26 = 1 boundary, if one makes reasonable allowance 
for the uncertainty in the Boron neutrino flux il- Finally, the global hts P, ^ have also 
found acceptable small angle MSW solutions for both these oscillation scenarios. But we do 
not consider them here, since the texture of our mass matrix naturally leads to bimaximal 
mixing as we shall see below. 


The dark matter constraint on the sum of neutrino masses comes from a recent global £t 
to the spectrum of density perturbation in the universe using various cosmological models 
T^ . The best fit is obtained with a hot and cold dark matter model, where the former 
constitutes 20 % of the critical density. Besides there is an astrophysical upper bound on the 
number of neutrino species from nucleosynthesis, which allows 1 or atmost 2 sterile neutrinos 
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We shall consider a four-neutrino mass matrix for the three doublet neutrinos and a 
singlet (sterile) neutrino. We shall present two scenarios, where the solar and atmospheric 
neutrino anomalies are explained by the oscillations (A) Vf. —»• vg and and (B) 

Vf. —>• Vr and > vg- The corresponding mass matrices will be related to one another 
by suitable permutation of neutrino indices. In each case maximal mixing between the 
oscillating neutrino pairs will be ensured by the texture of the mass matrix. Moreover we 
shall obtain the vacuum oscillation and the large angle MSW solutions in each case depending 
on the choice of parameters. 

(A) Vf. vg and v^ v^ Oscillations : 

In this case the texture of our neutrino mass matrix in the basis \Vf 




vs\ is 


m,f = 


/O 

0 

a 

\d 


0 

c 

h 

0 


a 

h 

0 

0 


d\ 

0 

0 

0 / 


( 1 ) 


Note that it has only 4 parameters. In comparison the earlier mass matrices considered 
had at least 5 parameters Moreover, the above mass matrix is minimal in 


the sense that it has only one diagonal element. The mass matrix of [^| has effectively 4 
parameters in the case of maximal vacuum oscillation solution of the solar neutrino. However 
it contains two equal diagonal elements, which could in general be different from one another. 
It is clear from the mass matrix that the neutrinoless double beta decay vanishes because 
Yfi Ulfrii = rriee = 0; so that the corresponding constraint |T^ is automatically satished. 

For the 3x3 submatrix of doublet neutrinos, we shall assume the hierarchy 


6 S> a, c. 


( 2 ) 


Consequently the and v^- will form a nearly degenerate pair with maximal mixing and small 
mass squared difference (~ 26c) to explain the atmospheric neutrino anomaly. Moreover, the 
remaining eigenvalue of this 3x3 submatrix gets a tiny double see-saw contribution 2a^c/6^, 
which will be much smaller than d over a wide range of the latter parameter. Consequently, 
the Vf. and vg will form a nearly degenerate pair with maximal mixing and small mass squared 
difference to explain the solar neutrino anomaly. The vacuum oscillation and the large angle 
MSW solutions will correspond to the choices d < a,c and d ~ 6 respectively. 
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I — Vacuum Oscillation Solution {b >> a,c > d) : In this approximation the mass 
eigenvalues are given by, 


mi 

m 2 

m 3 

m^ 


= d + 


a^c 


a?c 


c a? 

^ 2 ^ ^ ^ ^ ~ 262 
^ c a? 

^ 2~ W 


a^c 


= -d + 


a?c 


(3) 


and the corresponding mass eigenstates {yf = {z/i z /2 z /3 z/ 4 }) are related to the weak 

eigenstates (z/J = {ve t^s}) through the mixing matrix Uia as, 


( I2e\ 




Ur 


\l2s) 



Si 

J 


^ V2 


-Si 

1 

1 

V 2 

S 2 


51 
1 

1 

5 2 


Si 


\ 

1^2 
1^3 

T2 ^ 


(4) 


where, si = and we can neglect the terms S 2 and S 2 , which are of the order of ~ 
0(-^, |j) ~ 10“^ for our choice of parameters. For the given 4x4 mixing matrix Uia the 


' F262 ) 62 ) 

probability of two flavour oscillation is given by. 


2 (^rn^L' 


Puai^is ^a/s 4 ^ ^ UajUajUgjUgj sm . 

j>i \ 


(5) 


where, = mf — For our mixing matrix the flavour oscillation in each case is 


j' 


dominated by one mixing angle which can be determined by comparing the expression (^) 
with the effective 2 x 2 flavour oscillation formula 


Puai^/s sin 26ap sin 


'Am2.L' 

AE 


( 6 ) 


Thus we get an expression for sin^ 2,9ajs in terms of parameters of the mixing matrix Uia ■ 

For illustration we shall now present the solution for a specihc set of parameters, he., 
a = O.OSeV, b = l.SeV, c = O.OOleV and d = O.OOOleV. There are two pairs of nearly 
degenerate eigenvalues 


m^^ ~ —mt,^ ~ d = O.OOOleV 


~ —mj/g ~ b = l.SeV. 
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The LSND experiment can be explained by the oscillations between states with mass squared 
difference of the order eV^ which means that it can be explained by the oscillations between 
the 1 / 1^4 and 1 ^ 2,3 states. To explain LSND as an oscillation between the z/g and the effective 
mixing angle sin^ 206^ is obtained by comparing (^) and and reading off the mixing matrix 
elements from (^), 


sin^29e^, = -4:{UelUe2U^i U^2 U^\U-\- UQ 2 U^2 ^e4^e3^^4^/i3} 

- -4 X 4(-.;)(i)- 


( 8 ) 


Similarly the other masses and the relevant mass squared differences and the corresponding 
mixing angles for the experiments listed in Table 1 are given by, 




2a^cd 

= 


IP 


sin^ 26'es 

= 1 

= 

ml^ - ml^ 

= 2hc = 


sin^ 20^^ 

= 1 

^^\sND — 

- ^2 

= }p — d 


sin^ 29e^, 

= 8si = 

^DM = 

= 3eD. 


= 2.2 X IQ-^^eV 


4a^ 

Ip 


= 0.004 


( 9 ) 


The z/e —z /5 oscillation gives the vacuum oscillation solution to the solar neutrino anomaly, 
while the z/^ —oscillation explains the atmospheric neutrino anomaly. The LSND result 
is explained by z/^ —z/g oscillation. The contribution to dark matter is 3 eV. The CHOOZ 
0 and other laboratory constraints are satisfied. 


Note that the above solution consists of two nearly degenerate pairs of maximally mixed 
neutrinos, separated by a realatively large mass squared gap. This is known to be the 
favoured mass conhguration for satisfying the various laboratory constraints |^, The 
four model parameters are used to ensure that the three mass squared gaps correspond to 
the required values of Am^ for the solar, atmospheric and LSND neutrino oscillations, and 
Oefj, corresponds to the required mixing angle for LSND. The hot dark matter prediction 
comes out as a bonus. All these features are natural predictions of our mass matrix; and as 
such they will be shared by each of the alternative solutions discussed below. It should also 
be noted that the underlying double see-saw mechanism is responsible for generating mass 
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squared gaps differing by 10 orders of magnitude starting with mass parameters, which differ 
by only 3-4 orders of magnitude {i.e., similar to the case of the up type quark mass matrix). 

II — Large Angle MSW Solution {b > d » a,c) : In this approximation (6 ^ (i)the 
mass eigenstates are given by, 


mi 

m2 

m 3 

m4 


= d 


a?d 


+ 


a%'^c 


2(62 - d?) 2(62 _ ^ 2 ) 


'■b^c 


c c a 0 a' 

^ 2 ^ ^ ^ 2(62 - d 2 ) “ 2(62 - d 2)2 

c a% a%‘^c 

^ 2 ~ m ~ 2(62 _ ^ 2 ) “ 2(62 - d 2)2 


= -d + 


a?d a?b‘^c 


+ 


2(62-d2) 2(62-d2)2 


( 10 ) 


and the mixing matrix has the same form as (^, with si = j.^ 2 _)_^ 2 )i /2 > and S 2 = s '2 = ^^ 2 ) 1/2 

. Hence like before > vs mixing and the Vr mixing are maximal, where as the 

1^11 mixing is given by the small parameter si. It may be added here that one gets a 
smooth numerical solution at 6 = d, although the approximation ( 10 ) breaks down there. 

In this case let us consider a choice, a = 0.025eV, 6 = l.SeV, c = O.OOlSeH and d = 
1.25eH. Then the different masses and the relevant mass squared differences and the corre¬ 
sponding mixing angles are given by. 


m, 2 i ~ -mi .4 

~ d = 1.25el7 

my^ ~ -mi .3 

~ 6 = 1.5el7 

A 2 2 2 

Am,„; = 

2a%‘^cd 5 2 

= ( 4 . _ 

sin^ 29es 

= 1 

A 2 2 2 

= 26c = 0.0046^2 

sin^ 29 

= 1 

^'rnlsND = 

= 6 ^ - = 0.69el/2 

sin^ 29ef, 


^DM = 'E\^i\ 

= 5.5eV. 


i 


The numerical values of the mass square differences have been calculated using the exact 
solutions of for the masses. The analytical expressions for the mass square differences are 
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from the polynomial approximation (|^). The nnmerical agreement for the mass sqnare 
differences between the exact solntions and the polynomial approximation agrees upto 6 
decimals in this range of parameters. 

Thus the v^. —> vg oscillation provides the large angle MSW solution to the solar neutrino 
problem. The and z/g oscillations explain the atmospheric neutrino anomaly 

and the LSND result respectively as before. The contribution to dark matter is 5 eV. It may 
be added here that with these parameters Si < S 2 ; so the effective mass of z/g is slightly lower 
than that of z/ 5 . 

(B) z/g — v.^ and —> vg Oscillations : 

In this case we can use the same mass matrix as before if we make the following change 
of basis, 

(z/g Us) — ^ (z/g z/^ Us Ur). (12) 


Note that with this change of basis the diagonal element for the sterile neutrino mss con¬ 
tinues to remain zero, which is an important condition as we shall see below. Moreover the 
change of basis does not affect the mass eigenvalues mi, m 2 , m 3 and m^. But now the nearly 
degenerate pair mi and m 4 represent the two maximally mixed states of z/g and z/,-, while 
m 2 and m 3 represent similar admixtures of u^ and z/ 5 . Thus the solutions (^) and o will 
continue to hold with the exchange of the neutrino flavour indices r and S in des and d^r- 
Consequently they will represent the vacuum oscillation and large angle MSW solutions to 
the solar neutrino anomaly via z/g —Ur oscillation, while the atmospheric neutrino anomaly 
is explained via u^ —> us oscillation. The rest of the results remain the same as before. 

Let us briefly discuss the possible mechanisms underlying the above mass matrix. Con¬ 
sider hrst the 3x3 submatrix corresponding to the three left-handed doublet neutrinos. 
This sub-eV scale mass matrix could arise from the standard see-saw mechanism with three 
heavy right-handed singlet neutrinos. Alternatively, one can get it without any right-handed 


neutrino but with an expanded higgs sector via a radiative mechanism or Majorana cou¬ 
pling of the left-handed neutrino pairs to a heavy Higgs triplet p^. In both cases one can 
naturally obtain a sub-eV scale mass matrix. The extension of the mass matrix to include a 
light singlet neutrino has been tried recently in each of the above three models i 0 ,i- In the 
standard see-saw model it is assumed to be one of the right-handed singlets while one adds 
a singlet neutrino in the other two models. In each case one has to impose a zero Majorana 
mass for this singlet, as otherwise it will naturally assume a high mass value. This is the 
reason why we have set mss to zero in our mass matrix (1). In the standard see-saw model 










this has been done by assuming a singular Majorana mass matrix for the singlet neutrinos, 
so that one of the eigenvalues (mss) is zero p[. In the other two models the mss is made to 
vanish by imposing an additional symmetry p, |^. Finally one asks if this singlet neutrino 
can naturally have Dirac masses in the < 1 eV scale? It seems possible to get it in the Zee 
model 1^ and the triplet higgs model via the same suppression mechanisms which keep 
the 3x3 doublet mass matrix in the sub-eV range. But in the case of the standard see-saw 
model it had to be put in by hand |p. We feel it is important to look for a more natural 
way to keep the Dirac masses small in this model. 

To summarize, we have presented a texture of four-neutrino mass matrix which auto¬ 
matically ensures bi-maximal mixing between i^e ^ and —>■ or vice versa. Thus 
with only four parameters it can account for the solar, atmospheric and LSND neutrino 
anomalies while remaining consistent with other experimental constraints. The prediction 
of the desired hot dark matter density comes out as bonus. Depending on the choice of 
parameters we can get both the vacuum oscillation and the large angle MSW solutions to 
the solar neutrino anomaly. Thanks to the underlying double see-saw mechanism, one can 
generate the desired mass squared gaps differing by 10 orders of magnitude starting with the 
four mass parameters which differ by only 3-4 orders of magnitude. 
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